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Abstract—The hydrodynamic impact problem is posed in 
the frame of the linearized Wagner theory. The so-called flat 
disk approximation is used. The penetrating body is an elliptic 
paraboloid. The pure vertical entry is well known and exact 
solutions exist. When, after hitting the fluid, the body moves 
along its 5 other degrees of freedom, exact solutions are still 
available provided the angular motions are small enough. In case 
of a pure translational motions (both vertical and horizontal) it 
is possible to determine occurence of low pressure on the wetted 
surface either in free drop tests or constant velocity penetration. 


I, INTRODUCTION 


The linearized three-dimensional Wagner problem is known 
as a mixed Boundary Value Problem posed in Potential Theory. 
The linearization follows from the asumption that the relative 
body shape and free surface shape are flat enough. That is 
measured by the deadrise angle which cannot be smaller than 
4° and greater than 20°. Mathematically the Boundary Value 
Problem is posed in a half space, say z < 0, and the plane 
z = 0 is described by the cartesian coordinates (x, y), and 
it contains the free surface FS(¢) and the expanding wetted 
surface D(t), as illustrated in (Fig. 1) 


II, THEORETICAL DEVELOPMENTS, 


The flow is described with a velocity potential » which 
verifies the Laplace equation. The impermability condition on 
the wetted surface leads to a Neumann condition; the normal 
velocities of both the fluid and the body are continuous. On the 
free surface, by neglecting gravity and surface tension effects, 
the linearized Bernoulli equation reduces to a homogeneous 
Dirichlet condition. The Boundary Value Problem hence reads 


Ag=0 


z<0 
0, («.y) € FS(t) 
. (x.y) € Dit) 

(a? +1? +2?) > 00, 





(CD) 





The regularity of the solution is increased by using the 
displacement potential ġ which is the time integral of the 
velocity potential . The Boundary Value Problem written 
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Fig. 1. Illustration of the linearized three-dimensional Wagner problem. 


for the displacement potential remains unchanged except the 
Neumann condition on D(t), which now reads 


t 


b2 = —h(t) + f(x,y), nty= [ W(r)dr 2) 


o 
where the right hand side is the distance between a point on 
the body and the undisturbed free surface, that is to say its 
vertical displacement. Let us consider an elliptic paraboloid 
described in a Cartesian local coordinate system (O, X, Y, Z) 
as illustrated in (Fig. 2). Its equation is 

en Aa A 

Z= (KY) = o> +97 68) 





where (rz, ry) are the curvature radii. Alternatively we can use 
the coefficients A? = 2r, and B? = 2r,. We consider that 
the body penetrates the free surface not only with a vertical 
velocity but also with horizontal velocities along its two axis 
0 and Y 0 as described in (Fig. 3) We denote 
W) the Cartesian components of the body velocity. The 
motion is described in a coordinate system attached to earth 
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Fig. 2. Cut of the elliptic paraboloid by the vertical planes X = 0 (left) 
and X = 0 (right). 








Fig. 3. Kinematics of the elliptic paraboloid with vertical and horizontal 
velocity components. 


(A,a,y,2) which coincides with the Cartesian coordinate 
system (O, X,Y,Z) at initial time ¢ = 0. In the system 
attached to earth, the coordinates of a point on the wetted 
surface, are thus 


t 
r=X +f U(r)dr = X + dx(t) (4) 
0 


‘ 
y= +f V(r)dr = Y + dy(t), (5) 
o 


t 
z=Z =f W(r)dr = Z — h(t) (6) 


At any instant the vertical displacement of a point on the 
wetted surface is hence 





@—ôz(t))? | (y~dy(t))? 


=-h 
fa (O+ 2 


(7) 





v 
That means that the right hand side of the Neumann condi- 
tion remains unchanged provided that we use the change of 
variables 





x—dx(t), g=y—dy(t), ž=2z (8) 


The Boundary Value Problem is re-written for 4(#, 7, 
o(x.y.2,t) as follows 














Ad=0 z<0 

o=0 F=0, (7, EFS) (gy 
$z = h(t) + f(a) ==0, (#7) € Dit) 

670 (+9 +2) o0, 


and a standard boundary value problem is obtained. As proved 
in [3], the displacement potential is also the solution of the 
integral equation 


1 A2d dso 


x =— = h(t) — f(%,i), (10) 
27 Jot) ||M Moll 
with the additional boundary conditions 
ġ=0, and 6,=0, at T(t) = D(t)N FS(t). (11) 


The integrand of (10) is studied by [8]. In particular it is proved 
that Ago is singular as the inverse square root of the distance 
to P(t). In other words Ag@ can be written 


(12) 





where s is regular on D(¢) and the function F vanishes along 
the contact line I(t). It is conjectural whether or not F is 
the equation of an ellipse in the plane (%, 7) knowing that the 
shape is an elliptic paraboloid. Indeed, that result is proved in 
[6]: an elliptic paraboloid entering fluid with arbitrary vertical 
velocity yields an elliptic wetted surface as well. The function 
F can hence be expressed as 

zo 2 p 
F@i)=1-3-B (13) 
Then we can invoke Galin’s theorem (see appendix A) The 
application of this theorem to the present case is detailed in 


[7]. In particular a closed form expression of the displacement 
potential is obtained 

__2ha(, # p 2/8 

-3E a B 


That result proves that the displacement potential varies lin- 
early with the penetration depth. No special assumptions are 
made regarding the kinematics of the entering body. The 
expression of the velocity potential follows from 





(14) 


do 


=x 


which is consistent with the result of [4]. Galin’s theorem 
allows to prove that the aspect ratio of the elliptic wetted 
surface (noted k) is uniquely defined by the aspect ratio of 
the elliptic paraboloid (noted k-,), according to 


_# pit D/E 
Br, | 2-D/E 


whose variation is plotted in figure (Fig. 4). 





(15) 


k? (16) 
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Fig. 4. Variation of ky with k; see equation (16), 


Here Æ and K are Elliptic Functions of the eccentricity 





n/2 
o= f VI sin? 640 (7) 
A 
fe do 
Ko= f -== D=(K-E)/è (18) 
Kle) o yi- e?sin?0 ‘ ue: 
b=B,fh (e2) =ni (19) 
eT n(2- 2) (20) 


It is worth noting that the lengths of the semi-axes, at any 
instant of penetration, are proportional to the square root of 
the penetration depth h. Finally we can retrieve the expression 
of the displacement potential in the coordinate system attached 
to earth 


(x,y, t) = 





2ha (, (z-ôz)? (u - ôv)? \ 
Ze (1 EE (21) 


That result is essential in order to properly determine the 
pressure, 


III, HYDRODYNAMIC LOADS 
A, Pressure formulation 


The pressure follows from the linearized Bernoulli equation 


P= -ppa = -phe (22) 


The time derivatives in (22) make the expression of the 
pressure rather complicated. In order to illustrate its variation 
in space and time, we consider an elliptic paraboloid which 
enters the free surface with vertical velocity W = Å and the 
horizontal components U = dx and V = dy. Here are some 
details of the time derivatives 


VF (2 ` 
b= s (GTF + ha) 








(23) 


öz)? (y~6y)? 





F= -3 7 (24) 
a Uæ—ôr)U , Ay—dy)V 
Po + at 
2a(x— dx)? — 2b(y — dy)? 
4 2ale— day" is y) 25) 
By using the identities 
2ha = ha, £ = ; (26) 
it remains VF 
MVF 
=- 7 
— - ôy) V 
M=Wa+ aon gee + 2an t25 (28) 
To go further we need AY, yielding 
<a u? v? 
NI = aW — 2ah (= + 7) + 
-őz „ 2 os 
en (nù - Zy) + av) + 
Ay — » 2b = 
+205) (uw = Fv) + av) (29) 
and finally the pressure reads 
p= TF (“r + smë) (30) 


We can already anticipate that the pressure distribution over 
the wetted surface will be substantially different than for the 
pure vertical entry. 


B. Spatial and temporal variations of the pressure 

We focus on the variation of the pressure in space and 
time. More precisely we can determine the location of possible 
change of sign of the pressure. As an illustration, the quantity 
P = REYE is expressed along the axis of the y-translational 
motion. We reduce the motion of the entering body to sway 
U = 0. The coordinate Y = #620 is such that |Y] < 1 on 
the wetted surface. The pressure is expressed as follows 


P(Y) = Q(t, Y)(1 - Y?) + Xwv + bY) (Wa + 2hkYV) 





81) 
with 
j2 oy f 5 
ouv- wa- aea Bi (vasa -t 
(32) 
It is noticeable that at the contact point Y = —1 (ie. at the 





stern), the pressure vanishes if the velocity of expansion of 
the wetted surface b is set to the translational velocity V. That 
happens at a given instant. By using the expression (26) and 
assuming that V and W are constant, this instant is 


2w 2y 
n? W Bw D 
teeayao = Gyr = Gyr (+p 





(33) 
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with the notation of equation (19). The pressure can remain 
negative over a significant region of the wetted surface. As 
an example, we consider the elliptic paraboloid entering with 
constant velocities (V, W). We introduce m = = = ž 
which_appears as a nondimensional measure of time m = 
2 


7 





Th. The pressure (31) reduces to 
P(Y) =kVW (2mv2 +2+}- m) (G4) 
m 


(Fig. 5) shows the variation of {A with m and Y The 






ie __# 
Fig. 5, Variation of gbqy with m and Y according to equation (34) when 
velocities (V, W’) are constant in time, 


superposed line separates the plane (Y,m) in two areas, one 
of them corresponds to negative pressure. Above a critical 
coefficient merit, = 1/V2, P(Y) has zeroes and becomes 
negative. Figure (6) shows (in logscale) the variation of poqr 
(from equation 34) as m < Merit.. 





a 


' a 


Fig. 6. Logscale variation of giy with m < merit. and Y according to 
equation (34) when velocities (V, W) are constant in time. 


It is shown that the first zero appears at Y = —1/V2. Note 
that the Wagner model of water impact can be used only during 
the initial stage when the pressure is either positive everywhere 
in the contact region or negative strictly inside this region. 
After the time instant (given by equation 33)), when the zone 
of negative pressures approaches the contact line, Wagner's 
model cannot be used and another model should be employed. 
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C. Analysis of the pressure peak 


The pressure in (30) is said of first order. Its singularity 
behaves like the inverse square root of the distance to the 
contact line, where F = 0. The matching of this solution 
to a local solution at the spray root is commonly used to 
propose a uniformly bounded distribution of pressure. Then 
it is possible to capture its maximum. The twodimensional 
jet root solution is described by [13]. The asymptotics of the 
velocity potential (@,o.¢) and the pressure (Poor) along the 
contact line I require the introduction of velocity V; of fluid 
in the jet and the thickness of the jet d 


dE 





root = —AV5 Proot = 2pV} (35) 





7 


7 


where € is the distance between a point (a cosa, bsina) at 
the contact line and a point (XQ, Yq) on the wetted surface 


along the normal Ñ to the contact line 
I le ik 
IIVF\| Vecos? a + k?sin? a 


Figure (Fig. 7) illustrates those notations. As a first step, we 


cosa 


N=- A 
ksina 


ro 


a bo Y 


Fig. 7. Notations for the description of the flow at the contact line. 


use a Taylor development to express F in (23), F in (25) and 
M in (27) asymptotically at point Q close to the contact line, 
ie. as Ẹ is small. 


Masymp = Wa + 2h(U cosa + kV sina) +O(€) (37) 
Fasymp = Zy cos?a + k?sin?a +0(€°) (38) 

à w 2 P 
Pamp = q + ZU cosa +kV sina) +O(€) G9) 


Introducing these last expressions in (30) and (23) and match- 
ing the resulting velocity potential and pressure with (35) 
yields the fluid velocity in the jet Vj and the thickness of 
the jet d as follows 


1 Wa +2h(U cosa + kV sina) 


2h — Veos?a+k?sin? a 


Vy = (40) 
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2 


cos? a + k? sin? a)?/? 


(41) 
It is shown that the jet thickness does not depend on the 
horizontal motion. This matching also provides the measure 
of the pressure peak 


p (Wa+2h(U cosa +kV sina))? 


Pmax = 373 (42) 


cos? a + k? sin” a 


As a example, we determine this peak for translational motion 
along the Y axis (U = 0), yielding 
ʻe) 
eee 
m 


It is shown that the maximum is higher at the stern than at 
the bow part of the elliptic paraboloid. 


1 
pao = 40V? (i (43) 


D. MLM for the pressure 


MLM stands for Modified Logvinovich Model. It has been 
introduced by [9] in order to improve the prediction of the 
hydrodynamic force. However the equations can be rapidly 
heavy if we treat the general 6dof motion. That is why we 
keep horizontal translational motion only. 

As a starting point we express the impermeability condition 
on the body surface in terms of the velocity potential y which 
is the 3D solution describing the flow about the actual shape 


Voisin (44) 


where i is the exact normal unit vector on the wetted surface 
and 7 is the velocity of the body 


=-Wr+UE+Vi7 (45) 
Knowing the instantaneous body position, 
= S(w — dx, y — ôy) — h(t) = f*(z,y,t)—- h(t) (46) 


we hence determine the normal velocity on the body surface 
1 


2 (Pfa t Pufu Pa) 
yl+V7s 


Jp ñ= 


= Poa (Meme VS - Pa) (47) 
Vl+ VPP 
and the vertical velocity follows. 
2 =Ve- VE-t-Vf-W (48) 
We introduce the velocity potential ® 
D(x, y, t) = p(z, y, f* (x,y, t) — h(t), t) (49) 
from which we calculate the time derivative 
D= pe-p (ë-7s+W) (50) 
and its spatial derivatives 
F = ampt pS (51) 


yielding 


p= 5 (70 Ws—a-Vs-W), 


We get the final expression of the reduced pressure 
P(x, y,t) = p(a,y, f*(x, y,t) — h(t), t) 


D=1+V"F (52) 


—1 P(e,y,t) = O,+ lg- la +Ü? (53) 
P 2 2 , 


It remains to calculate the solution ® from the Taylor devel- 
opment 
ylz. y, 2, t) = p(z, y,0,t) + 29,:(£,y,0,t) + 
at z= f*(x,y,t)— h(t) 


where p(x, y,0,t) is the solution of the linearized Wagner 
model noted ip and its vertical derivative on the body surface 
¢.2(2,4,0,t) = p(x, y,0,t) = -W 


(x,y, t) =p - w[i t) -a)] 


(54) 


(55) 


where y") is given by equation (27). Its time and spatial 
derivative are then expressed as 





=. MF 
M= s HE) (56) 
MVF 
Vo) = VFYM + 57 
T WF 6n 
which are used in 
de= pP -W -h) +W FSW) (58) 
o = Ve —WVF (59) 


mae the MLM pressure with horizontal velocity reads 


-2 P(e, wt) =e? —W(st —h) + iw? + 39 


-3 (We -J vs)’ 


In order to compare the pressure (60) to the first order pressure 
(30), we consider a practical case of an elliptic paraboloid 
with curvature radii rą = VA = 0.75m and ry = VB = 2m 
respectively in the two horizontal directions x and y. It enters 
vertically with velocity W = 1m/s and the horizontal velocity 
along the y axis varies from V = 0m/s up to V = 10m/s. 
The pressure is measured with the sensor of finite size. Its area 
is circular with radius rs = 0.0025m. Figures (Fig. 8) and 
(Fig. 9) show the influence of the horizontal velocity on the 
pressure at two points along the axis of motion, |Y | = 0.05m 
(Fig. 8) and |Y| = 0.20m (Fig. 9). The effect of the pressure 
averaging over the sensor area is clear. Without this averaging, 
we would get rather unphysical sharp peaks. It is shown that 
the pressure is higher on the upstream side than along the 
downstream side. The case with a zero horizontal velocity 
V = 0m/s is intermediate between the pressures measured at 
both sides. It is worth reminding that the pressure maximum 


(60) 
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Fig. 8. Time variation of the pressure measured by a pressure sensor located 
at [Y | = 0.05m along the major semi-axis. The sensor is circular with radius 
ra = 0,0025m, 
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Fig. 9, Sce caption in (Fig. 9); the pressure sensor located at |¥| = 0.20m 


is evaluated without reference to the inner solution as in (42). 
The fact that dynamics of expansion is not affected by the 
horizontal motion, appears clearly since the duration of the 
pressure rise-up, as it is not instantaneous, is similar whatever 
the considered side and whatever the horizontal velocity. 

The distribution over the wetted surface is computed ac- 
cording to MLM. The influence of the horizontal velocity is 
illustrated in the figures (Fig. 10) to (Fig. 12) for V = 5m/s, 
V = 10m/s and V = 20m/s. We distinguish positive and 
negative pressure in the plane (m, Y) with different colors. In 
addition we multiply the pressure with VI — Y? so that the 
change of sign of the pressure appears clearly. For a given size 
of the elliptic paraboloid, that is to say before separation of the 
flow occurs, the smaller the horizontal velocity, the latter the 
negative pressure appears. The superimposed line corresponds 
to the theoretical limit given by zeroes of the first order 
pressure (34). The discrepancies between first order pressure 
and MLM pressure seems to be greater at low horizontal 
velocity, 


E. Force 


At the leading order, the vertical force F; follows from 
the integration of the pressure along 7. We can check that 
performing the integration of the pressure over the contact 

















iation of the MLM pressure in terms of m and Y’ for a horizontal 
5m /s. The superimposed line: zero of the first order pressure 








Fig. 11, See caption in (Fig, 10), Horizontal velocity V = 10mn/s, 


region D yields the same result that 


1? 
viJ fst 


since both ġ and its normal derivative on the contour &D(t) 

vanish. By using the expression of the displacement potentiel 

in (21) and the identities (26), the vertical force (whatever the 
4np a Are 2,3 d? 


entering kinematics) is 
hanana = 15/2 
I5E d? TE" " ae [ ] 


It is worth noting that the force depends on time variation of 
the penetration depth h only since the factor © ‘a constant 
and depends on the geometry only. It is rather surprising that 
the force does not depend on the pressure distribution while 
the latter strongly depends on the horizontal velocities and 
varies significantly during the penetration. 

Assuming that the vertical velocity W is constant, the force 


reduces to 1 
Broke we Vi 


Figure (Fig. 13) compares experimental time variation of the 
force data (see [10]) and the formula (63) scaled by W°/? As 
usual the first order formulation (63) overestimates the force. 


(61) 


[abn] (62) 





F, = (63) 
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Fig. 12, 








Fig. 13. Time variation of the force sealed with W°/2. Comparison with 
experimental data deesribed in [10]. 


When the kinematics is a pure vertical motion the computation 
of the force from a numerical integration of the pressure 
MLM is detailed in [9]. The corresponding force is plotted 
in (Fig. 13) as well and a better agreement with experimental 
data is obtained, 

It is worth reminding that the numerical integration is 
better handled if the position of the vanishing pressure is 
precisely known in the close vicinity of the contact line. In 
the present case, it is shown that the area where the pressure 
is negative starts growing after a while as the body penetrates 
the fluid. No attemps are made to compare the first order force 
(62) to the numerical integration of the pressure MLM when 
the translational velocity exists. However in the experiments 
mentionned above, the penetration before separation of the 
flow is far smaller that the time at which b = V. In that 
case MLM results with or without horizontal velocity are very 
similar, 

It is also worth noting that the experimental vertical force 
starts increasing with a finite slope at initial instant of contact. 
It is suggested to average the force over a short time duration 
thus accounting for the fact that force sensor cannot measure 
instantaneously. 

Considering now a free drop of the body. Its mass in air is 
H. In the absence of any other hydrodynamic force that the 


force (62), Newton's law yields a differential equation which 
can be integrated. The following equation must be solved 


ph + C(e)h*/? = pWot (64) 
where W, is the initial velocity and the factor C(e) 
- e)n? 
Cle) = U e 65) 


15E(e) 


is a constant only function of the geometry in terms of the 
aspect ratio k according to (16). In order to illustrate the 
kinematics during free drop, it is easier to plot the time ¢ in 
terms of the penetration depth h(t) for given data (jz, e, Wp). 
The following figure shows such variation 





Wont merona Apana soana 






The added mass of the penetrating body reads 


5 ` 
Ha = gel? (66) 
and the acceleration reads 
z 15 vh 
À = -=C le) W 
h Tole WWor a (67) 


From equation (30), it is possible to define other circumstances 
when the pressure calculated at the initial contact point X 
Y =0, may vanish during free drops, That occurs when AV 
0, yielding 





U? + K?V? = aa W (68) 

By using Newton's law (uW = —Wrya), and the identity 

jia = E after some algebra, it is shown that if the 

horizontal kinematics is fitted to the vertical kinematics as 

follows 

W2a? p — 2jia 
2h p+pa 








U? + RPV? = (69) 
the pressure at the initial contact point vanishes thus leading 
to possible cavitation. If the horizontal velocity is zero, the 
vanishing pressure occurs when the added mass is twice the 
mass in air of the body. That result is shown in [11]. Further 
investigation of a growing cavitation pocket is described in 


12). 
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APPENDIX A 
GALIN’S THEOREM 
Theorem: /et the functions (G, F, P, I) defined as follows: 


1 1 z y? 





Oleo) = i= aF F(z,u)=1- 37 
(70) 
M M-p 
Pleu) =} J vpz”! m) 
p=0 420 
P(x, 
Ilev) (x0; yo) 


G(x — 10, — yvo)drodyo (72) 
p VF (xo, yo) i Z 


then, if (x,y) € D={(x,y), +B <1), 
M M-p 


Ie) =Y Y apr”! 


=0 4=0 


(73) 





that means that I is a polynomial in (x,y) of the same degree 
as P. 

We assume here that Galin’s theorem can be inverted: Zf I is 
polynomial in (x,y), say with degree M, then P is polynomial 
as well with same degree M. The expressions of coefficients 
Ayq in terms Of bpg are detailed in [1]. 
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